Cooper Pair Shape in Normal-metal/Superconductor Junctions 
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In s-wave superconductors the Cooper pair wave function is isotropic in momentum space. This 
property may also be expected for Cooper pairs entering a normal metal from a superconductor 
due to the proximity effect. We show, however, that such a deduction is incorrect and the pairing 
function in a normal metal is surprisingly anisotropic because of quasiparticle interference. We 
calculate angle resolved quasiparticle density of states in NS bilayers which reflects such anisotropic 
shape of the pairing function. We also propose a magneto-tunneling spectroscopy experiment which 
could confirm our predictions. 
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It is well known that Cooper pairs consisting of two 
electrons are characterized by electric charge 2e, macro- 
scopic phase, internal spin, and by time and orbital struc- 
tures m . The charge 2e manifests itself in various exper- 
iments, like Shapiro steps, flux quantization and excess 
current due to the Andreev reflection. The macroscopic 
phase generates the Josephson current [l|. The inter- 
nal spin structure is classified into spin-triplet and spin- 
singlet states which can be identified by nuclear magnetic 
resonance. Further, based on a symmetry with respect 
to the internal time, superconducting state can belong 
to the even-frequency or the odd-frequency symmetry 
class i3|] • The orbital degree of freedom is described by an 
angular momentum quantum number I. In s-wave super- 
conductors with I — 0, the Cooper pair wave function is 
spherically symmetric on the Fermi surface, i.e. an angu- 
lar structure of a Cooper pair is isotropic in momentum 
space. The shape of Cooper pairs in p-wave (/ = 1) and 
d-wave {I = 2) superconductors is characterized by two- 
fold and four- fold symmetries, respectively [^. The well 
established properties listed above hold in bulk supercon- 
ductors. The presence of perturbations like spin-fiip or 
interface scattering may change the symmetry of Cooper 
pairs. For instance, unusual odd-frequency property of 
Cooper pairs inproximity structures was predicted in re- 
cent studies [J, y . The shape of Cooper pair wave func- 
tion in non-uniform systems like superconducting junc- 
tions is not necessarily the same as that in the bulk state. 
Despite the extensive study of the proximity effect dur- 
ing several past decades, rather little attention has been 
paid to the problem of Cocroerpair shape in non-uniform 
superconducting systems [y, |7| . This issue is quite im- 
portant in view of current interest to the physics of su- 
perconducting nanostructures. 

The aim of the present Letter is to clarify the con- 
sequences of breakdown of translational symmetry in 
superconductors on the Cooper pair shape. For this 
purpose, we study the proximity effect in quasi two- 
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FIG. 1: A schematic illustration of a N/S junction. 



dimensional normal metal / superconductor (N/S) junc- 
tions by solving the Eilenberger equation, treating self- 
consistently the spatial variation of the superconducting 
pair potential. We analyze the pairing function and the 
local density of states (LDOS) in N/S junctions with 
spin-singlet s-wave and spin-triplet p^^-wave supercon- 
ductors. Surprisingly, the pairing function in a normal 
metal turns out to be strongly anisotropic even in junc- 
tions with s-wave superconductors. To detect the com- 
plex Cooper pair shape, we propose to use scanning tun- 
neling spectroscopy in rotating magnetic field. We show 
that the calculated tunneling conductance exhibits com- 
plex patterns even in the s-wave case. 

Let us consider a quasi-two dimensional N/S junction 
as shown in Fig. [1] which is the simplest example of non- 
uniform superconducting system, where the S region is 
semi- infinite and the normal metal has finite length L. In 
the case of pa;-wave superconductor, we assume for sim- 
plicity that triplet Cooper pairs consist of two electrons 
with opposite spin projections on z-axis, i.e. Sz — 0. 
These assumptions do not limit the generality of the dis- 
cussion below. We consider a perfect N/S interface with 
full transmissivity, while it can be shown that charac- 
teristic behavior of Cooper pairs remains qualitatively 
unchanged even in the presence of a potential barrier at 
the N/S interface. 

The quasiclassical Green's functions [81] in a normal 
metal (N) and a superconductor (S) are parameterized 



5± — /i±^i + J2±^2 +g± T-3, i3± j — -L, Uj 

where a superscript i(= iV, 5) refers to N and S, fj 
(j = 1 — 3) are the PauH matrices, and 1 is a unit ma- 
trix. The subscript +(— ) denotes a moving direction of a 
quasiparticle in the x direction [8|, and A+(x) (A_(a;)) is 
the pair potential for a left (right) going quasiparticle. In 
a normal metal, A±(x) is set to zero because the pairing 
interaction is absent there. The Green's functions can be 
expressed in terms of the Ricatti parameters |9| , 

/« = TMrT{^) + d\^Mi + r'-^{^)C^±\^)i (2) 
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-1 for i ~ N . The param- 
eters T\^ (x) and Ct' {x) obey the Eilenberger equation of 



with Vi ~ \ for i ^ S and m; 



the Ricatti type [9|, 

«wf.a,r«(a;) = -A±(a;)[l + (r«(a;))2] + 2eiy,r^i\x), 
iVF,d^C^^\x) = -A±(a;)[l + id\x)f] 



2eiy,C^^\x), 



with £ = £ -|- i(5o where e is the energy of a quasiparticle 
measured from the Fermi level and Sq is the inverse of the 
mean free time due to impurity scattering. In the clean 
limit, we consider 6o <?C Aq. Boundary condition at a; = 



.(N) 



(N), 



~L is given by Cj_ (~L) = — F-^ {—L). Boundary con- 
dition at the N/S interface becomes C±^^(0) = -F^^^(O) 

and C±^^(0) =_-F^^^(0). The pair potential A±(a;) is 
expressed by A±(x) = A(a;)$±(6')9(x), where a form 
factor $±(0) is given by <&±(6') = 1 for s-wave symme- 
try and ± cos for p^^-wave one with 6 being an incident 
angle of a quasiparticle measured from the x direction. 
Bulk pair potential is A(oo) — Aq, and we determine the 
spatial dependence A(a;) in a self-consistent way. 

For a; ^ io, the angular structure of /jjl follows 
that of the pair potential, whereas f^^ is zero with 
Lq = 'VfI'^'^Tc being a coherence length and Tc being 
the transition temperature. The pairing function j^^ is 
generated by inhomogeneity in a system and thus has 
a finite value only near the interface and in a normal 
metal. Recent study [lOJ showed that /}2 has an odd- 
frequency symmetry. Generally speaking, functions /{^ 

and J2-L have opposite parities. If a superconductor has 
s- (Px)-wave symmetry, the induced odd- frequency com- 
ponent has the odd (even) parity, respectively. Pair- 
ing function /i is defined in the angular domain of 
-7r/2 <Q< 37r/2. We denote /i(6') by /i+(6l) in the 
angle range -7r/2 <Q < it 12 and fi{9) = /i-(7r - 9) for 
7r/2 < < 3n/2. The angular structure of functions /2 
and g is defined in the same manner. LDOS is given by 
the relation pL{d) = Keal[g{d)]. In what follows, we fix 



temperature T = O.OSTc, the length of the normal region 
L = 5Lo, and So = O.OIAq. 

In Figs. [2] and [3] we show polar plots of /i, /2, and 
PL in a s-wave and p^-wave junctions for several choices 
of £ and X. Black, red and blue lines represent, respec- 
tively, the results for a; = oo (superconductor), a; = 
(interface), and x = —L/2 (normal metal). The odd- 
frequency component is always absent for a; = oo. Since 
PL is independent of x in N, the resulting value of pL at 
X = is equal to that at a; = —L/2. 

First, we focus on the s-wave case (Fig. [5]). At £ = 
in Fig. [2ja), the even- frequency component /2 has a 
circular shape reflecting s-wave symmetry for x = oo 
and 0. However, at x = —L/2, the shape is no longer a 
simple circular but a double distorted circles. The shape 
of /2 in superconductor always has the circular shape 
independent of £ as shown with black broken lines in 
Fig.2(a), (d), and (g). At e = O.IAq in Fig.2(d), /s at 
the interface (red line) slightly deviates from the circular 
shape, while the shape in N drastically changes. The 
tendency is more remarkable at £ = 0.5Ao in Fig. 2(g). 
The butterfly-like pattern seen in /2 at the interface (red 
line in Fig. 2(g)) is completely different from the original 
circular shape in a superconductor. At £ = 0, /i at the 
interface becomes ellipsoidal as shown in a red line in 
Fig. 2(b). The shape of /i at x = and x = —L/2 shows 
the butterfly-like pattern (Fig. 2(e)). For e = 0.5Ao, 
the line shape of /i has many spikes as shown in Fig. 
2(h). Such anisotropic property of /i and /2 affects the 
LDOS as shown in Fig. 2(c), (f), and (i). In particular, 
LDOS at the interface for e — O.SAq (red line in (i)) 
strongly deviates from the circular shape. The LDOS 
in a superconductor vanishes at £ < Aq (black lines in 
Fig. 2(f) and (i)). At the interface, the shape of pL in 
Fig. 2(c) is quite similar to that of /i in Fig. 2(b). These 
profiles can be qualitatively understood as follows. At 
X = 0, the following relations hold 
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eyi-p[2ie L/{vf cos 6)]. For s <C Aq, F ^ 1/i and fi± ~ 
ig± are satisfied. Thus shape of function /i is similar 
to that of PL- This argument seems to be valid even 
for £ — 0.1 Aq in Fig. 2(e) and (f), and for e — O.SAq in 
Fig. 2(h) and (i). The oscillating behavior in /i, /2, and 
PL is more remarkable at £ = 0.5Ao. Although we do not 
present calculated results of /i and /2 at x = —L/2 for 
£ = 0.5Ao, the butterfly-like pattern with many spikes in 
the pairing functions can be seen also in a normal metal. 
The directions of the spin projections in LDOS are char- 
acterized by small value of S, which has close relation to 
the formation of the bound states 1]|. For ~ ±7r/2, 
a oscillates rapidly with small variation of 6, which ex- 
plains the fine structures in LDOS around = ±7r/2. 
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FIG. 2: (Color online) The results for s-wave symmetry. The 
shape of the even- frequency pair amplitude [(a),(d),(g)], the 
odd- frequency one [(b),(e),(h)], and the angle resolved local 
density of state [(c),(f),(i)]. Red lines: s = (at N/S inter- 
face), blue lines: x = —L/2 (in a normal metal) and black 
lines: a; = oo (in a superconductor), e = for (a), (b) and 
(c), e = 0.1 Ao for (d), (e) and (f), and e = 0.5Ao for (g), (h) 
and (i). The angle 6 is measured from the x-axis. 



The quasiparticle interference effect is a source of bound 
state formation in a normal metal. As a result, the cir- 
cular shape of Cooper pairs in s-wave superconductor is 
modified into the butterfly-like pattern in a normal metal. 
Next, we discuss the results for pa;-wave junctions 
shown in Fig. [3l fn a superconduct or (a: = cc), f unc- 
tions /2 and g are given by Aq cos 9/ sj e^ — Aq cos^ d and 
ej \Je'^ — Ag cos^ 0, respectively. As shown by black lines 
in Fig. 3(d), (f), (g) and (i), the amplitudes of fi and g 
become large along the directions Q — cos~^(e/Ao). At 
e = and x = 0, formation of a mid-gap Andreev reso- 
nant state [13, ll3| significantly enhances the amplitudes 
of /i and ph compared to that of /2 as shown by red 
lines in Fig. 3(a), (b), and (c). For e = and e — O.IAq, 
the shapes of /i and pL at the N/S interface are similar 
to those in s-wave superconductor junctions (red lines in 
Fig. 3(b), (c), (e) and (f)]. The shape of /i and J2 in the 
N region is rather complex too. At e = 0.5Ao, /i, /2 
and ph also exhibit the butterfly-like patterns as shown 
by red lines in Fig. 3(g), (h), and (i). Similar to s-wave 



FIG. 3: (Color online) The results for p^-wave symmetry. The 
notations are the same as in Fig. [2] 



case, functions /i and ]i in N have a complex line shapes 
with many spikes. The butterfly-like pattern in angle re- 
solved local density of states could be observed through 
the LDOS in magnetic fields. 

Here, we propose an experimental setup of scanning 
tunneling spectroscopy (STS) in the presence of the mag- 
netic field. As shown in Fig. [I] magnetic field is ap- 
plied parallel to the N/S plane. Tunneling current at a 
fixed bias voltage is measured as a function of the an- 
gle (/) between the x axis and the direction of magnetic 
field. The vector potential in this configuration is given 
by {A^.Ay) 



14, 151. We 



-XH exp(— z/A)(sin 0, cos </)) 
assume that thickness of a quasi two-dimensional su- 
perconductor is sufficiently small compared to the mag- 
netic field penetration depth A. Magnetic field shifts the 
quasiparticle energy e to e — -ff Aq sin((/) — 0)/ Bq where 
Bq = /i/(2e7r^^A) and ^ — fivp /■k/S.q. For typical values 
of ^ ~ A ~ lOOnm, the magnitude of Bq is of the or- 
der of 0.02Tesla. Local density of state observed in STS 
experiments is given by, p{(j)) = J_^i2PL{9,<p)d9. 

In Fig. m p is plotted as a function of </>. For e = 0, 
AqH/Bq gives the effective energy of a quasiparticle. In 
an s-wave superconductor, it is evident that p is indepen- 
dent of (j) and the amplitude of LDOS is vanishingly small, 
see curve A in Fig. HJa). In a p^-wave superconductor, 
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FIG. 4: (Color online) Magneto-tunneling conductance as a 
function of the orientation angle of magnetic field in a spin- 
singlet s-wave junction [(a) and (c)] and in a spin-triplet px- 
wave one [(b) and (d)]. e = in the cases (a) (b) and e = 
0.5Ao in the cases (c), (d). A: a:: = cx) (in a superconductor) 
for H ^ O.lBo, B: X ^ (at N/S interface) for H = 0.05Bo, 
and C: x^Oior H = O.lBo- 



as shown by curve A in Fig. [D^b), p depends slightly on 
(f> due to quasiparticle excitations near the nodal points 
at the Fermi surface [l6|. At the N/S interface, p in the 
s-wave case has minima at = tt/2 and 3tt/2 and max- 
ima at ^ = and tt, as shown by curve B in Fig. |4ja). 
On the other hand, in the pa;-wave case, p has maxima 
at (/) = 7r/2 and 37r/2, see curve B in Fig. \M^>). This 
behavior can be explained in the following way. Eq. ^ 
yields approximate expression for p at H ^ Bq 



tt/2 



-7r/2 



2(1 -a2) de 
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i^ ± 2a cos[C sin (6* — (f))/ cos ( 
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where a = eyip{—AL6o/h\vFx\), C = 4:LH/h\vFx\, and 
the sign in the denominator is + (— ) for s-wave (px-w&ve) 
junctions. We note that a is a positive number almost 
independent of 9 and C is a small positive number. For 
(j) — riTT + 7r/2 with integer n, the magnitude of the ar- 
gument of cosine function in the denominator becomes 
small and the denominator is reduced to (1 ± a)^ for s- 
wave (j)a;-wave) case. When (f> deviates from mr + n/2, 
the cosine function decreases and therefore p has a dip 
(peak) at = tt/2 and 37r/2 for s-wave (p^-wave) cases 
as shown by curve B in Fig.|3I^a) ((b)). Similar argument 
also explains the maximum (minimum) oi p at (j) — m: 
in s-wave (p^^-wave) junctions. When the field H is in- 
creases further (curves C in Figs. 4(a) and (b)), the shape 
of the p exhibits more complex behavior reflecting the 
butterfly-like pattern of pL shown in Figs. [5] and [31 be- 
cause the increase of H has qualitatively similar effect as 
the increase of e. 

The complicated oscillating features are seen much 



more clear at e = O.SAq. As shown by curves C in Figs. 
4(c) and (d), the period of oscillations becomes shorter 
than that in Figs.HJa) and (b), reflecting the butterfly- 
like patterns in Figs. [2ji) and Figs. Elji). At the same 
time, the magnitude of the oscillations at £ = O.SAq be- 
comes smaller than that at e = 0, since the integration 
with respect to 9 averages the butterfly-like pattern in 
PL- The discussed features in LDOS at the N/S interface 
differ strongly from those in a bulk superconductor. We 
conclude that STS experiments in magnetic field should 
resolve the remarkable deformation of Cooper pairs. 

In summary, we have studied the Cooper pair shape 
in normal-mctal/supcrconductor (N/S) junctions by us- 
ing the quasiclassical Green's function formalism. The 
quasiparticle interference leads to striking deformations 
in the shape of a Cooper pair wave function in a normal 
metal. As a consequence, the angle resolved local den- 
sity of states exhibits the butterfly-like patterns. We also 
show that the anisotropic shape of Cooper pairs could be 
resolved by scanning tunneling spectroscopy experiments 
in magnetic field. The Cooper pair deformation is a com- 
mon feature of non-uniform superconducting systems in 
the clean limit. This provides a key concept to explore 
new quantum interference phenomena in superconduct- 
ing nanostructures. 
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